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Abstract 

jrt ' We propose a diffusion model for the recently discovered diffusion-induced Ramsey narrowing 
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arising when atoms diffuse in a buffer-gas cell in the laser radiation field. The diffusion equation 
for the coherence of metastable states coupled with an excited state by laser radiation of different 
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^ ■ 

Q,' frequencies in a three-level scheme of the atom-field interaction is obtained in the strong-collision 



approximation. The dependence of the shape of an absorption line near the transmission maximum 
of one of the frequencies on the two-photon resonance detuning for various geometries of the cell 
is investigated. 
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I. INTRODUCTION 

In a three-level system subjected to two laser fields coupling two metastable states (or a 
Hietastable state and a stable one) with an excited one, a dark or light-nonabsorbing state, 
namely a coherent superposition of two metastable states, can be formed. This phenomenon 
is called a coherent population trapping (CPT). The condition of the formation of a CPT 
state is the two-photon resonance under interaction of an atom with light, where the dif- 
ference between the frequencies of two laser fields is equal to the frequency of a transition 
between metastable states. If this condition is realized, one observes an abrupt decrease of 



M- 



Due to the CPT 



the fluorescence intensity of the atom in the laser radiation field 
phenomenon, it is possible to create a window in the absorption spectrum. As a result, 
light can propagate almost without losses through a medium that absorbs light under usual 
conditions, which represents the well-known phenomenon of electromagnetically inducec 
transparency (EIT) [4]. In addition, dark resonances are used for the light slowing down |5| 
and in the construction of compact laser frequency standards 6|] and underlie the effective 
method of population transfer between different states of atoms or molecules - stimulated 
Raman adiabatic passage (STIRAP) [7(]. The resonance width depends on both the coher- 
ence decrease rate for the lower states and on other factors, particularly on the pattern of 
atomic motion in a buffer-gas cell. The latter aspect will be the focus of our attention in 
this work. 

In the case where atoms move through a laser beam of finite width, the role of the 
coherence time is played by the residence time of an atom in the field. If the cell contains 
a buffer gas in addition to active atoms, their residence time in the laser beam increases. 
As a result, narrow resonances with a width of the order of tens of hertzs are registered in 
buffe.gas cells g The autho. of wo.ks gQ e„,pha..e the .o.e of the buffe. gas in the 



experiments on coherent population trapping in a three-level system, though it is considered 
that, after atoms have left the region of interaction with radiation, they do not return there 
anymore. 

A more detailed description of the process of atom-field interaction must take into account 
that, having left the region of interaction with the field, the atom can return there again 
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12|, so that atoms can interact with radiation several times before the loss of coherence. 



Thus, diffusion of atoms in a buffer gas essentially affects their response to the resonance 



FIG. 1: Three-level system: il.a and ilh are the Rabi frequencies 



excitation by the laser field. If the coherence relaxation time of the metastable states 
considerably exceeds the time till the repeated atom-field interaction, the atom can get back 
without loss of coherence after having spent some time in a dark region (beyond the beam). 
As a result, it is worth expecting the narrowing of the resonance line. This phenomenon 
w. called a diffu=ion-Wuced R.„.ey na..owmg Q (by analogy with the Itat,.ey method 
of separated oscillating fields [IJ, |15|). 

The diffusion- induced Ramsey narrowing of the transmission spectrum in the case of EIT 
observations was investigated in [l3| adducing the results of experiments with rubidium 
vapor and neon serving as a buffer gas, as well as the results of theoretical calculations. 
According to [l3|, the necessity of taking the diffusion process into account depends on the 
laser-beam diameter - the contribution made into the signal by atoms that had spent some 
time beyond the beam and got back becomes determinative with decrease in the beam di- 
ameter. The same authors subsequently published a detailed description of the developed 
theory 16|. Experimental results and theoretical calculations unambiguously confirm the 
physical interpretation of the phenomenon of diffusion- induced Ramsey narrowing of the 
transmission spectrum observed in [l3|. In particular, with increase in the laser-beam diam- 
eter or the buffer gas pressure, the shape of the spectral line changes from the non-Lorentz 
to Lorentz one in accordance with a decrease of the contribution made by atoms that have 
come back from the region beyond the laser field. In the cited works, the motion of atoms is 



described in the form of the Ramsey sequences: in each of them, an atom spends the time 
t']^ in the radiation field, moves in the dark region during the time ^2^*, returns back to the 
radiation region for the time t™, and so on. In order to find the transmission spectrum, the 
density matrix describing the Ramsey sequence was integrated with the probabihty distribu- 
tion, i.e. one performed averaging over all possible trajectories. The probability distribution 
was obtained in 16[] from the diffusion equation. 

We propose an alternative approach to the description of the diffusion-induced Ramsey 
narrowing. It is based on the diffusion equation obtained from the initial motion equations 



for the density matrix on the basis of the strong-collision approximation 



171]. The trans- 



mission spectrum will be directly obtained in the equilibrium approach instead of averaging 
over trajectories. As compared to works [13] and [l6|, we consider the relaxation at cell 
walls (supposing that the coherence is broken due to a collision of an atom with the wall) 
and take a realistic (Gaussian) intensity distribution in the radial direction. In addition, 
we also take the finiteness of the gas cell in the direction parallel to the laser beam into 
account. In order to understand how the dimension of the problem influences the result, we 
also consider the one-dimensional case where the cell is infinite in the direction of the laser 
beam, while the beam itself is infinite along one of the transverse coordinates. 

II. BASIC EQUATIONS 

Let us consider a gas of three-level atoms with the excited state ]2) and the metastable 
lower states ]1) and ]3). A field with the frequency Ua couples the states ]1) and ]2), while 
that with the frequency Ub couples the states ]3) and ]2) (see Fig. 1). The interaction of these 
fields with atoms is described by the Rabi frequencies Qa = M12 ' ^a/^ and Qi, = ^1^2 ' ^t/^, 
respectively. As the beams a and b are spatially bounded, the Rabi frequencies depend on 
the position r of an atom in space. The wave vectors will be considered close in magnitude: 

iU r^^ h>j r-^ hi 

l^a — i^b — i^- 

The equations for nondiagonal elements of the density matrix in the rotating-wave ap- 
proximation have the form 

iQ* 

P12 = i{Aa - kv)pi2 + -^ (P22 - Pll)- 



P23 = -^(Afe - kv)p23 + -^Pl3 + 



i^b , s , fdp22,\ , /5P23 



P31 = -z(Aa - Af,)p3i - — P23 + ~^12 + 

Here, the terms with the index "coll" describe the relaxation processes due to collisions of 
active atoms (interacting with the field) with atoms of the buffer gas resulting in a change of 
the velocity, while the coherence is conserved. The terms with the index "relax" describe the 
rest of relaxation processes, A^ = UJ21 — oJa and A^ = uj2z — cob are the one-photon detunings, 
and Aw = Aq — A^ is the two-photon detuning. It is assumed that 

dpi2\ _ -p fdp2A _ 

J- 12P12, ^, — —J- 23P23, 



•^^ / relax V "^^ / relax 

where ri2 and r23 stand for the coherence relaxation rates for the transitions 2—7-1 and 
2 — )■ 3, respectively. If Fis and r23 are sufficiently large, the collision terms in Eqs. (P) 
and (|2]) can be neglected. The relaxation rate for the forbidden transition |1) — )■ |3) is 
non-zero due to collisions between atoms, 



dp3i 



•^^ ' relax 



-713P31- (4) 



In the general case, the expression for the collision term 

'dpij{^r,v',t)\ , 



18|-|20| can be presented as 



dt 



+ y"ir,,(v',v)p,,(r,v',t)rfv', (5) 

where i^jj(v',v) is the collision kernel, and z/ is, in the general case, a complex- valued 
quantity with the frequency dimension. It can be interpreted as the collision frequency in 
the case where the scattering amplitudes in both states i and j are identical (see 18|). 



The collision term can be sim plifi ed using the strong-collision approximation (light atoms 
are scattered by heavy particles 17|]) and assuming that v and Kij are real, while the kernel 
i^jj(v', v) does not depend on v', i.e. the velocity of an atom v after a collision does not 
depend on its velocity v' before it. In this case, the velocity distribution (arbitrary in the 
general case) turns into the Maxwellian distribution after only several collisions, that is, an 

n 

atom quickly forgets its initial velocity. As was shown in [ITj], i^jj(v) = vW{y\ where Vr(v) 
is the Maxwellian distribution. Thus [20|, expression ([5]) takes the form 

'dpi 



dt 



•[p,^-W{v)N,,] 



coll 



where Nij{r,t) = / dv pij{r,v,t) 



(6) 



can be interpreted as the number of atoms with values of pij lying in the unit volume in the 
neighborhood of the point r, and u denotes the collision frequency 18|. 

In this approximation, collisions with the buffer gas change only the external degrees of 
freedom of atoms. 

Using the approximation |^| ^ ri2, r23 and neglecting the collision terms in ([1]) and (|2]), 
one can find the stationary solutions for pi2 and p23- After that, with regard for the fact 
that the nondiagonal elements pu are close to the equilibrium values pli (they differ from 
the latter by a small quantity of the second order in the field intensity), one obtains the 
equation for pai: 



P31 



-(7i3+iAw)p: 



K^a 



'31- 



+ 



ntn 



b"b 



T23-i{Ah-kv) Ti2 + i{Aa-kv) 



Psi 



"a"fe 



(0) 
Pll 



(0) 
P22 



+ 



(0) 

P33 



(0) 

P22 



Ti2 + i{Aa- kv) T23-i{Ab~kv) 



dp. 



31 



dt 



(7) 



coll 



The second term on the right-hand side describes the field broadening. It can be neglected 
in the case of weak fields and large relaxation rates Tij . 

In order to simplify the equation, we also neglect the Doppler broadening in the third 
term considering that T^j ^ Aq^^, Au, kv. 



Taking into account that the equihbrium elements of the density matrix pa are propor- 
tional to the distribution function W{\), we obtain a kinetic equation of the Boltzmann 
type for pgi. 

In the stationary case of interest, it has the following form: 

(v ■ V)p(r, v) = -(z/ + 7 + iAc<;)p(r, v) 

+W{w)[X{v) + uN{v)\. (8) 

Here and below, we use the notations 

p(r,v) = p3i(r,v), 7 = 7i3, 



2/„2 
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where N^i is determined by ([6]), and Wq is the normalization constant of the Maxwellian 
distribution, 

Xir)^\[^^^^P^y.ni. (9) 

Here, we took into account that P22 ^ Pn , P22 ^ P33 • ^^ °^^ "^^^ ^^^ from (|9]), the 
function A(r) ~ i7a(r)n^(r) describes the transverse profile of the beams. 
For the sake of simplicity, we introduce the substitution 

tto = 7 + i^^ , a = u + 'J + iAu = u + a^ . 

It is assumed that the collisions of atoms with walls result in the failure of the coherence 
between the lower states, so that 
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P(r,v)|rgs = 0, 

where 5* denotes the surface confining the buffer-gas cell. 

The shape of the spectral line is determined by the function T(Au;) = Re [S'(A(X')/S'(0)], 
where S{Auj) has the form 

S{Auj) = ff dr dw A(r)p(r, v) = /" A(r)A^(r) dv . (10) 

In order to find it, it is necessary to obtain the function p. 



In the case of a laser beam with the Gaussian intensity distribution in the plane normal 
to the direction of its propagation, the expression for A(r) in the cylindrical coordinates 
takes the form 

A(r,y.,^) = Aoe-'-'/'^% (11) 

where Aq is determined by i^, and r lies at the beam axis. 

In the strong-collision approximation, the time between collisions Tj^ = 1/z/ is small as 
compared to the characteristic time of flight of an atom through the interaction region r^ = 
a/vQ. Therefore, we consider that UTa ^ 1. 

In what follows, we investigate the effect of the beam size and the distance to the walls of 
the gas cell on the lineshape specified by the function ReS'(Aa;). The form of this function 
depending on the dimension of the problem will be considered as well. 

III. SOLUTION FOR INFINITE REGION 

For an infinite region, Eq. (|8]) can be solved with the help of the Fourier transformation 
with respect to the argument r: 

p(k, v) = /" p(r, v) e-'^' dr , iV(k) = f p(k, v) rfv . 

Hence, 

where the function 



F(k) = / ^^^ = # e-V^^^g erfc f -f- 
J a + lii.'v \k\vo \\k\vo 



\k\voJ 

describes the Voigt profile. Let us find the signal 5*00 using the properties of the Fourier 
transformation: 

SU^uj) =JdrN{r)X{r) = ^ J dkN{k)X{k) , (13) 

where n denotes the space dimension of the region (1, 2, or 3). Substituting A^ from (IT2ll . 
we obtain 

^ ' (27r)"y l-z/F(k) 



In some cases, the value of S{Au) is mainly determined by small k due to the factor X^{k). 
That is why the function F{k) can be replaced by the asymptotic expansion in the quadratic 



approximation (see 



): 



It is valid for sufficiently small k, for which |A;|fo/z/ ^ 1. 

A. One-dimensional case 

In the given case, A(x) = Aqc"^ /° , and ( !T^ has the form 



2, 



2 J ao + u^/td 

— oo 

2 ao Vy2; ' ^ ^ 

where 

2 2a;ao 2aa? 2 2 

(the quantity r^j will be interpreted below). In our approximation, v ^ 7, Aw, that is why 
one can consider T£, = const. The right-hand side of (IT5l) can be presented in the form of a 
superposition of the profiles 

00 
Sg\Auj)^ [ s{Auj,u)gg\u)du, (16) 



where 

s{Au,u} = . . .. , 7eff(M) = 7 + M /rz) 

7cff(M) +«Au; 

with the weighting factor (independent of u) 

£H-) = -4;;Te""^/^ (17) 

The combination r('u) = T£,/u^ can be called the effective diffusion time of an atom; t{u) 
for "fast" atoms is larger than that for "slow" ones. 



The profile Sio (Aco) can be described as the effective Lorentzian with the center at the 
origin of coordinates 

Sl{Aco) = ^° ■ 
i + iAui 

Its width To is determined by the relation 

_25(0) 

° S"{0)' 

Simple calculations for •yTu <^ 1 yield 

r^ = ^7^-fi-/fv7^+...V (18) 

Therefore, we obtain Fq — a/|7 for small 7. 



B. Two-dimensional case 

Using the similar procedure for the function A(r) = Aqc"^ /" in the polar coordinates, 
we obtain 

TcXW f ue~^ /^ du 



S^^\A 



OJ 



^0^ 



2 7 ao + u^/td 



^ e^^aV2 Ei, (^ ) . (19) 



Here, Eii denotes the integral first-order exponent 23|. The obtained expression is similar 



to flT5|) and has the same interpretation, though with the weighting function 

g^^{u) = ^ue--''\ (20) 

The width of the effective Lorentzian for '-jtb ^ 1 has the form 

r^ = 7^2 [In 2 - 7s - ln(7rB)] + o{^td) . (21) 

IV. EFFECTIVE DIFFUSION EQUATION 

A. Derivation of diffusion equation 

According to f llOp . the complex signal S'(Aw) can be expressed with the help of the zero- 
order momentum A^(r) of the distribution p(r) with respect to v. Let us find the equation 
for the moments of higher orders N'^'^\v). 
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First, we consider the one-dimensional case. Here, 

oo 



Equation ([8]) presented in the one- dimensional case as 

^M^ = -ap{x, v) + W{v) [X{x) + z/AtW (x)] (22) 

will be multiplied by v'' and integrated over v. This procedure yields a chain of equations 
for the moments N^''\x). 

Writing down the obtained equations separately for even and odd k and successively 
substituting the following equation into the previous one m times, we obtain the following 
expression for A^*^"): 

™ ' ,2A:\ ^2fc 



'"• I 2k\ J2k 



fc=0 

]_ ^2m+2jY(2m+2) 
Q,2m+1 ^rj,2m+2 ' ^ '' 

The result for higher space dimensions will be similar, though more complicated as the 
moments will be specified by tensors. A similar procedure yields 

A;=0 



,N^'\r) = Y, 



Oil , v^ , / „, 



x(V,...V.,J[A(r) + z/iVW(r)]+... (24) 

Expressions fl23|) and fl24l) represent asymptotic expansions, where one can leave only the 
first terms in the case of sufficiently large v (and a). In particular, the second-order terms 
result in the diffusion equation. 

Using the known expressions for the averages (viVj), one obtains 

\ na^ J 
in the quadratic approximation, where n stands for the space dimension. In view of the 
equality (f^) = ^Vq, one derives 

«oiV(r) = ^ AN{r) + A(r) , N{oo) = (25) 
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for the arbitrary n (accurate to terms of the order of l/v). This equation can be interpreted 
as the diffusion equation with the absorption coefficient ao = 7 + ^Aw and the complex- 
valued diffusion coefficient 

n.^4. pa) 

The diffusion coefficient is identical for all space dimensions. 
Expression ( !26l) at large v becomes real and turns into 

D = ^. (27) 

Proceeding from the formula a = a/Utd, we obtain the characteristic diffusion time (ap- 
proximate time, for which an atom leaves the beam) 



^D = T^ = -^. (28) 

-^0 



D vl 



For the further consideration, it is convenient to put down the diffusion equation (125|) in 
the form 

AiV(r) - /32iV(r) = -/(r) , (29) 



where /(r) = ^A(r), /3^ = ^"°° ~ ^^- Solving it with the help of the Green function 
method, one can see that the expression for the signal obtained by solving ( 129|) in the case 
of large v will be the same as that obtained earlier by direct calculations f lT^ . One can 
also see that the number of terms used in the asymptotic expansion -F(k) of expression (TT^ 
correlates with the number of terms of Eq. (1241) that must be kept in order to obtain the 
diffusion equation. 

In the case of a finite cell, the kinetic equation ([8]) can be solved formally, by interpreting 
the last term on the right-hand side as a nonuniform one (see 22]). The result will be the 
same as that derived from the solution of (1251) accurate to the terms z/^^ (to which the 
diffusion equation is actually valid). 

The atomic motion is characterized by five characteristic times: 
_ a _ R _ 1 

''"a ) 'O? ) ''"7 ) 

vq vo 7 

1-^ = - , rD = — = iyT^. (30) 
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According to the accepted approximations, they satisfy the following conditions: 

Tu <Ta<TR, Td <T./, T^ <€. T^ , Ta < Td ■ (31) 

The first condition results from the strong-collision approximation and the geometric con- 
figuration R > a for small z/ and (or) large R. The approximation used for the derivation 
of the diffusion equation will be valid for a finite-size cell. The second and third conditions 
correspond to the slowness of relaxation processes, while the last condition is evident. 

Let us introduce the dimensionless time and space scales t = 7t and f = r/a and denote 
the dimensionless velocity by -Oo = ^0/(70), the diffusion coefficient by -D = D/{'~fa'^), the col- 
lision frequency by i> = Z//7, and the cell size hjR = R/a. The dimensionless characteristic 
times will have the form 

1 2i) 

I R 

ra = ira = —, TR=-iTR = —. (32) 

As will be seen from the further solution, the shape of the transmission line in the infinite 
case is completely determined by two dimensionless characteristic times f^, and fo, whereas, 
in the case of a finite cell, it depends on its size - R (width) and / (length, in the three- 
dimensional case). 

The diffusion equation in the form (!29l) will be solved for different space dimensions. 
The most effective approach is that based on the Green functions. In order to simplify the 
comparison of the cases of finite R and i? — > oo, we present the result of investigating the 
effect of space confinements on the signal in the form most similar to the solution for an 
infinite cell. 

B. One-dimensional case 

Using the general solution of the one- dimensional diffusion equation with the boundary 
conditions N{±R) = 0, we obtain 

+00 ^ 
^(1) ^^^ ^^ f Xiu/a)Xniu/a)du 



-00 
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where \R{u/a) = \{u/R)b^^\u, R), and b^^\u,R) is a factor depending on R: 
b^^\u,R) = B^^\u/a,R), 



Xlk) J 



R 

ch(/3x) 



m_^ 



—ikx —ikR 



ch{/3R) 



dx. (34) 



This expression means that, similar to the infinite case ( jTBj) . Sp^ can be written down as a 
superposition of Lorentzians: 



+00 
S^j^\Auj)r^ f s{Auj,u)g^j^\u,Auj)du, 



though, in contrast to flT7|) . the weighting function in the case of finite R depends on R 

and Au: 

g^^\u,Auj) = gg\u)b^'\u,Au,R); 

moreover, Mmji^oo 9^ (u, Alu) = goo{u). In the case of i? — )■ 00, the dependence on Au 
disappears. The factor g}^ {u,Au) cannot be interpreted as a weighing one due to the 
dependence on Au. In addition, the Lorentzian s{Au, u) cannot be replaced by the more 
complicated profile sr{Au, u) = s{Au, u) ■ b^^\u, Au, R) in order to separate out the weight, 
as it was done in the infinite case, because this "profile" sr becomes singular for some values 
of u (in particular, for u — ?■ 00), though this singularity is compensated by the other factor 
goo (u). Thus, the interpretation of S}^ (Au) as a weighted superposition is impossible here. 
For the Gaussian intensity distribution, the function b^^\k, R) can be expressed in terms 
of the error functions: 



,, . , erfc (^ + i^) + erfc (^ - i^) 
5(1'^) (A;, R) = l ^-^ ^-^ ^-^ ^ 



2 

(,.^,.)/, cos kR erf (f + f) + erfc (f - f ) 

^h^ 2 • ^^^^ 

In the case i? ^ a, -^ ^ ^, one can obtain the asymptotic behavior of B^^\k,R) for 

small k: 

5(1'^) {k, i?) ^ 1 - e'^'^'+^'^^'Z' ^^ . (36) 

chpir 
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The graphs of ReSu^Au) and the associated parameters for the typical values 

7 = 1.0x102Hz, z/ = 1.0 X lO^Hz, 

a = 1.0 X 10~^ m , vo = 2.95 x 10^ m/s, 
and td = 2.3 x 10^^ s with the corresponding dimensionless values 

td = 2.3 X 10"^ , % = 1.0 X 10"^ 

are given in Fig. 2. 

One can see that the line will be narrower for larger R (where the probability to return 
without loss of coherence is higher). Formally, the geometric factor gj^ {u, Aco) for finite R 
is larger than g^{u). Therefore, the profile ReSR{Au) will be wider. 

C. Two-dimensional case 

Using the axially symmetric solution of the two-dimensional diffusion equation in the 
region r < R with the boundary condition A^(r) | _ „ = 0, one obtains 



where Xr^u/o) = \(u/a)b^'^\u, R) and 
b^^\u,R) = B^^\u/a,R), 



^"'*-«>-|i/^^"- 





(2), 



R 



dr ■ (38) 



moreover, livsiR^oobR {u, R) = 1. 

The same way as in the one-dimensional case, the weighting factor of the Lorentzian has 
the form fif^j (m, Auj) = gio {u)b^'^\u, R, Au). The interpretation of fl38l) and the comparison 
with the infinite case (fT9l) are the same as those in the one-dimensional one (135|) . The graphs 

(2) 

of ReS')j (Acj) in Fig. 3 are plotted for the same values of the parameters as in Fig. 2. 
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^/a = infinity a 
^/a = 10 
^/a = 20 
^/a = 30 




Aco/y 




FIG. 2: One-dimensional case: normalized transmission T = Ke S (Aoj) / S (0) as a function of 
Aoj/j for different i? (a) and the transmission spectrum halfwidth as a function of R/a (b). The 
values of the parameters are given in text 

D. Three-dimensional case 



In the three-dimensional case where A^(r) is zero on the cylindrical surface of radius r = R 
confined by the planes z = ±/, the diffusion equation is identical to that obtained in the 
two-dimensional case, whereas the signal has the following form (modified two-dimensional 
profile) : 



S'n'i^^) 



IGvra^ ao 



■X 



X 



do + uVtd 



(39) 
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RIa = infinity 
f?/a = 1 
RIa = 20 
RIa = 30 




FIG. 3: Two-dimensional case: normalized transmission T = ReS'_R(Aw)/S'_R(0) as a function of 
Auj/j for R = oo and finite R (a) and the transmission spectrum halfwidth as a function of R/a 
(b). The values of the parameters are the same as in Fig. 2 

where \{k) denotes the two-dimensional Fourier transform of A(r), b^'^\u, R) is the factor 
from the two-dimensional case depending on R (see ( 155]) ). and 

tanh^(u, Au) 



b^^\u,l,Auj) = l 



au,Au) 



^{u,Aio) = -^/u^ + /3^a^. 



(40) 



The weighting factor of the Lorentzians in the three-dimensional case (!39|) has the form 

= g^^\u)b^'\u,R,Au)b^'\u,l,Aco). 
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i?/a = infinity, //fl = 20 
■ RIa = 10, lla = 20 
- RIa - infinity, lla - 50 

,R/fl=10, //a = 50 




Aco/y 



^(3), 



:(3)/ 



FIG. 4: Three-dimensional case. T = ReS'jj (Au;)/S']j (0) as a function of Aio/'y for different 
space parameters 




FIG. 5: Comparison of the functions T„ = Ke Soo (Auj) / Soo (0) for n = 1,2,3 



W, 



:-(")/ 



Its properties are similar to those in the one- and two-dimensional cases. It is worth noting 
that its dependence on the cell dimensions R and / appears in the form of independent 
factors. 

The transmission spectra normalized to one at the maximum for different values of the 
space parameters are given in Fig. 4. As expected, a decrease of / results in the broadening 
of the spectrum. Figures 5 and 6 present the transmission spectra for all dimensions in the 
case of an infinite cell and a cell with the size R/a = 10. One can see that the narrowest 
spectrum is obtained in the one-dimensional case and the widest - in the three-dimensional 
one. The values of the parameters used for plotting the graphs in Figs. 4-6 are the same as 
in Fig. 2. 
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FIG. 6: Functions T„ = Re s''^\Au}) / S^ (0) , n = 1,2,3, for R/a = 10 




FIG. 7: Normalized transmission T = Re Si, (Au) / Si, (0) — comparison with experiment 



V. COMPARISON WITH EXPERIMENT 



Figure 7 shows the comparison of the results of numerical calculations with experimental 
data [l3]. The calculations were performed for the dimensionless parameters 

f„ = 3.39 X 10"^ , fr = 1.69 X 10-3 , 



r^ = 1.00 X 10-' 
t)o = 2.95x 10^ 
^ = 5 , 1 = 5 



td = 2.30 X 10" 



D = 4.35 X 10' 



close to the experimental conditions. One can see that our theoretical results are in rather 
good agreement with those of experimental measurements. 
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It is worth noting that the experimental data presented in Fig. 7 also agree well with 
the theoretical calculation of the transmission spectrum obtained by averaging over atom 
trajectories. This testifies to the fact that the proposed model is also in good agreement 



with that based on the averaging over trajectories 13|, ll6|, at least in the cases where 
these theories can be consistently compared. The proposed model has a wider range of 
applications, since the model with averaging over trajectories can be used only for the case 
of an infinite cell (at least in its present form). 

VI. CONCLUSIONS 

We have constructed a model for the description of the phenomenon of Ramsey diffusion- 
induced narrowing of the transmission spectrum of atoms in a buffer-gas cell for the case of 
weak fields in the strong-collision approximation. This model can be used for an arbitrary 
intensity distribution of laser beams in the direction normal to their propagation. The 
general theory is illustrated by calculations for the Gaussian intensity distribution. 

The analytical expressions for transmission spectra obtained with the help of the effective 
diffusion equation qualitatively agree with experimental data and the results obtained by 
averaging over atom trajectories. The proposed model gives a possibility to investigate the 
shape of the transmission spectrum as a function of not only the intensity distribution in the 
plane of the laser beam but also the size of the buffer- gas cell. We have considered different 
geometric configurations (one-, two-, and three-dimensional). Comparing the spectra for 
one-, two-, and three-dimensional models, one can see that the line becomes wider with 
increase in the dimension. 

The work is carried out in the framework of Projects F28. 2/035 and RFFD/1-09-25. 
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